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The DuBois-Reymond Fundamental Lemma of the 
Fractional Calculus of Variations and an Euler-Lagrange 
Equation Involving only Derivatives of Caputo 
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Abstract 



Derivatives and integrals of non-integer order were introduced more than three centuries 
ago, but only recently gained more attention due to their application on nonlocal phenomena. 
In this context, the Caputo derivatives are the most popular approach to fractional calcu- 
lus among physicists, since differential equations involving Caputo derivatives require regular 
boundary conditions. Motivated by several applications in physics and other sciences, the 
fractional calculus of variations is currently in fast development. However, all current for- 
mulations for the fractional variational calculus fail to give an Euler-Lagrange equation with 
only Caputo derivatives. In this work, we propose a new approach to the fractional calculus 
of variations by generalizing the DuBois-Reymond lemma and showing how Euler-Lagrange 
equations involving only Caputo derivatives can be obtained. 
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1 Introduction 



The fractional calculus with derivatives and integrals of non-integer order started more than three 
centuries ago, with l'Hopital and Leibniz, when the derivative of order 1/2 was suggested |L. This 
subject was then considered by several mathematicians like Euler, Fourier, Liouville, Grunwald, 
Letnikov, Riemann, and many others up to nowadays. Although the fractional calculus is almost 
as old as the usual integer order calculus, only in the last three decades it has gained more attention 
due to its many applications in various fields of science, engineering, economics, biomechanics, etc. 
(see [2HS] for a review). 

Fractional derivatives are nonlocal operators and are historically applied in the study of non- 
local or time dependent processes. The first and well established application of fractional calculus 
in physics was in the framework of anomalous diffusion, which is related to features observed in 
many physical systems, e.g., in dispersive transport in amorphous semiconductor, liquid crystals, 
polymers, proteins, etc. [TIE] ■ Recently, the study of nonlocal quantum phenomena through frac- 
tional calculus began a fast development, where the nonlocal effects are due to either long-range 
interactions or time-dependent processes with many scales 14( 9-12]. Relativistic quantum me- 
chanics [T5Hl5j . field theories [TMT8] . and gravitation [TI5], have been considered in the context 
of fractional calculus. The physical and geometrical meaning of the fractional derivatives has 
been investigated by several authors [201EI] • It has been shown that the Stieltjes integral can be 
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interpreted as the real distance passed by a moving object [201111] • The most popular operators 
of fractional calculus, namely the Riemann-Liouville [T] and Caputo 22, 23 fractional operators, 
are then given similar physical interpretations |20[|21) . 

One of the most remarkable applications of fractional calculus in physics is in the context of 
classical mechanics. Riewe showed that a Lagrangian involving fractional time derivatives leads 
to an equation of motion with nonconservative forces, such as friction 24, 25 . It is a striking 
result since frictional and nonconservative forces are beyond the usual macroscopic variational 
treatment and, consequently, beyond the most advanced methods of classical mechanics |26j . Riewe 
generalized the usual variational calculus to Lagrangians depending on fractional derivatives, in 
order to deal with nonconservative forces [24II25] . Recently, several approaches have been developed 
in order to generalize the least action principle and the Euler-Lagrange equations to include 
fractional derivatives 24,25,27 34 . In this new formalism, because of the nonlocal properties 
of fractional time derivatives, the Euler-Lagrange equations appear to not respect the causality 
principle. This difficulty is under investigation by several authors [35]|36] . On the other hand, in all 
the above mentioned references, the Euler-Lagrange equations always depend on the Riemann- 
Liouville or mixed Caputo and Riemann-Liouville derivatives. For the state of the art of the 
fractional calculus of variations and respective fractional Euler-Lagrange equations, we refer the 
reader to the recent book [37] . 

It is important to remark that while the Riemann-Liouville fractional derivatives pQ are histo- 
rically the most studied approach to fractional calculus, the Caputo [2"2"ll2"3"] approach to fractional 
derivatives is the most popular among physicists and scientists, because the differential equations 
defined in terms of Caputo derivatives require regular initial and boundary conditions. Further- 
more, differential equations with Riemann-Liouville derivatives require non-standard fractional 
initial and boundary conditions that lead, in general, to singular solutions, thus limiting their 
application in physics and science [BJ. On the other hand, within current formulations of the 
fractional calculus of variations, even Lagrangians depending only on Caputo derivatives lead to 
Euler-Lagrange equations with Riemann-Liouville derivatives (see [38] and references therein). 
This is a consequence of the Lagrange method to extremize functionals: application of integration 
by parts for Caputo derivatives in the Gateaux derivative of the functional, relates Caputo to 
Riemann-Liouville derivatives. In the present work, we propose a different approach, by generali- 
zing the DuBois-Reymond fundamental lemma to variational functionals with Caputo derivatives. 
This new result enable us to obtain an Euler-Lagrange equation in integral form, containing only 
Caputo derivatives. Moreover, when the Lagrangian is C 2 in its domain, we obtain an Euler- 
Lagrange fractional differential equation depending only on Caputo derivatives by taking the 
Caputo derivative on both sides of the Euler-Lagrange equation in integral form. 

The article is organized as follows. In Section [3J we review the basic notions of Riemann- 
Liouville and Caputo fractional calculus, that are needed to formulate the fractional problem of 
the calculus of variations. The fractional DuBois-Reymond lemma is proved in Section |3l while 
in Section [4] we obtain the Euler-Lagrange equations in integral and differential forms. Finally, 
Section[5]presents the main conclusions of our work. Although we adopt here the Caputo fractional 
calculus, the results can be straightforward rewritten within other approaches, like the Riemann- 
Liouville one. 

2 The Fractional Calculus 

There are several definitions of fractional order derivatives. These definitions include the Riemann- 
Liouville, Caputo, Riesz, Weyl, and Grunwald-Letnikov operators [TH5l l39ll40] . In this Section, we 
review some definitions and properties of the Caputo and Riemann-Liouville fractional calculus. 
Although our first main objective is to obtain a DuBois-Reymond lemma for functionals with 
Caputo fractional derivatives, the Riemann-Liouville derivatives are also defined here because 
they naturally arise in the formulation. 

Despite many different approaches to fractional calculus, several known formulations are con- 
nected with the analytical continuation of the Cauchy formula for n-fold integration. 
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Theorem 2.1 (Cauchy formula). Let f : [a, b] — > R be Riemann integrable in [a, b]. The n-fold 
integration of f , n G N, is given by 



f{x)(dx) n = II I •••/ / f(xi)dxidx 2 ---dx n -idx n 

(I) 

1 r m du , u 



T(n)J a (x-w) 1 -" 

where T is the Euler gamma function. 

The proof of Cauchy's formula ([T]) can be found in several textbooks, for example, it can be 
found in p] . The analytical continuation of ([T]) gives a definition for integration of non integer (or 
fractional) order. This fractional order integration is the building block of the Riemann-Liouville 
and Caputo calculus, the two most popular formulations of fractional calculus, as well as several 
other approaches [TH51 l3T)PI)] . The fractional integrations obtained from (JT|) are historically called 
Riemann-Liouville fractional integrals. 

Definition 2.1 (Left and right Riemann-Liouville fractional integrals). Let a £ K+. The opera- 
tors a Jx an d xJb ' i defined on L\([a,b]) by 

•™ = fi)f(^ 4 < 2 > 

and 

where a, b £ M with a < b, are called the left and the right Riemann-Liouville fractional integrals 
of order a, respectively. 

For an integer a, the fractional Riemann-Liouville integrals ([2]) and ([3]) coincide with the 
usual integer order n-fold integration ([T]). Moreover, from definitions ([2]) and ([3j) , it is easy to see 
that the Riemann-Liouville fractional integrals converge for any integrable function / if a > 1. 
Furthermore, it is possible to prove the convergence of ([2]) and ([3]) even when < a < 1 [40] . 

The integration operators a J" and x J£ play a fundamental role in the definition of Caputo and 
Riemann-Liouville fractional derivatives. In order to define the Riemann-Liouville derivatives, we 
recall that, for positive integers n > to, it follows the identity D™f(x) = D^ a J^~ m f(x), where 
is the ordinary derivative d m /dx m of order m. 

Definition 2.2 (Left and right Riemann-Liouville fractional derivatives). The left and the right 
Riemann-Liouville fractional derivatives of order a G M + are defined, respectively, by 

a D a x f{x) := D n xa J n x ~ a f{x) 

and 

x D%f{x) := (-l) n D^ x ,rr a f{x) 

with n = [a] + I; that is, 

•ww = Tx^y t [ ( x-!)L-^ w 

n« fn (-1)" d n f b f(u) 

The Caputo fractional derivatives are defined in a similar way as the Riemann-Liouville ones, 
but exchanging the order between integration and differentiation. 



and 
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Definition 2.3 (Left and right Caputo fractional derivatives). The left and the right Caputo 
fractional derivatives of order a G M+ are defined, respectively, by G D f{x) := a J x ~ a D x f(x) and 
%D£f(x) := (-l); l J"- a L>J/(x) withn=[a] + l; that is, 



On«^._ 1 r / (n) («) 



z!D f(x) := — r / . ' (6) 

a x,/w T{n-a)J a (x-m) 1+q -" w 



and 

Dhf(x) := 7T7 — ■ — r / t — - — rrr" (7) 

6 w r(n - a) J x (u - x) 1 + a ~ n K ' 

where a < x <b and f( n \u) = d £ Li([a,b\) is the ordinary derivative of integer order n. 

An important consequence of definitions ([!"]) — (ffj) is that the Riemann-Liouville and Caputo 
fractional derivatives are nonlocal operators. The left (right) differ-integration operator Q and 
(|6]) (([!} and ©) depend on the values of the function at left (right) of x, i.e., a < u < x (x < u < b). 

Remark 2.1. When a is an integer, the Riemann-Liouville fractional derivatives (|4]) and ([5]) 
reduce to ordinary derivatives of order a. On the other hand, in that case, the Caputo derivatives 
© and differ from integer order ones by a polynomial of order a — 1 [31140] . 

Remark 2.2. For < a < I, and because they are given by a first order derivative of a frac- 
tional integral, the Riemann-Liouville fractional derivatives (j4]) and ([5]) can be applied to non- 
differentiable functions. One can even take Riemann-Liouville derivatives of nowhere differ entiable 
functions, like the Weierstrass function |41j . In fact, one of the main reasons for using fractional 
calculus is the possibility to study non- differentiability [421143] . 



Remark 2.3. R is important to note that, while the Caputo derivatives ^ and ^ of a constant 
are zero for all a > 0, „ D 1 = ~. D b 1 = 0, the Riemann-Liouville derivatives ((4]) and ([5]) are not 
zero for a £ N: 

(x~ar a _ (b-x)-<* 

r(i-a) ' x h r(i-o) 

For aGN, the right hand sides of (JSJ) become equal to zero due to the poles of the gamma function. 
Furthermore, for the power function (x — a)* 3 one has 



.^i = ^zr. ^l = ^l_. (8) 



a DZ{x-af = v J^ l) A x-aY-«. 
T(p + 1 - a) 

Remark 2.4. Let a G]0, 1[. If the Riemann-Liouville and the Caputo fractional derivatives of 
order a exist, then they are connected with each other by the following relations: 

G a D%f{x) = a D«f(x) - =£^(x - a)-° 
1(1 — a) 

and 

c x DU{x) = x D a b f{x) - ^f^—Ab x)- a . 

T(l - a) 

Therefore, if f(a) = (f(b) = 0), then c a D a J(x) = a D a J(x) (%D°f(x) = x D%f(x)) [HQ]. 

We make use of the following three properties of the Riemann-Liouville and Caputo derivatives 
(Theorems 12 .21 [2~3l and [2~4]) in the proofs of our results. 

Theorem 2.2 (see, e.g., [3"|l40]). Let a > 0. For every f E Li([a,b]) one has 

a LY* xa rj{x) = f(x), x D« x J a b f{x) = f{x) (9) 



and 

almost everywhere. 



^Dl a rj{x)=f{x), %D b!t Jtf(x) = f(x), (10) 
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Theorem 2.3 (Fundamental theorem of Caputo calculus — see, e.g., [3,40]). Let < a < 1 and 

f be a differentiable function in [a, b] . The following two equalities hold: 

a J% C a D a J{x) = f{b)-f{a) (11) 

and 

b J a a CD a b f(x) = f(a)~f(b). (12) 

Theorem 2.4 (Integration by parts — see, e.g., [39]). Let < a < 1 and f be a differentiable 
function in [a, b]. For any function g 6 Li([a, 6]) one /las 

" g(x)%D«f(x)dx = [ f(x) x D%g(x)dx + [ x D<£- l g{x)f{x)\ \ (13) 



b g(x)ZD%f(x)dx= f f{x) a D^g{x)dx+[{-l) n a Dt- 1 g{x)f{x)] b a . (14) 

The proof of relations ©-([121) can be found in several textbooks — see, e.g., [31140] . Relation 
^ is a direct consequence of Q and |JSJ|. Theorem 12.31 is the generalization of the fundamental 
theorem of calculus to the Caputo fractional calculus. It is important to mention that pip and 
(|12[) do not hold in the Riemann-Liouville approach. Finally, we remark that the formulas of 
integration by parts (Tj"3"]) and (fT4"|) relate Caputo left (right) derivatives to Riemann-Liouville 
right (left) derivatives. 



3 The Fractional DuBois— Reymond Lemma 



The fundamental problem of the fractional calculus of variations consists to find a function y that 
maximizes or minimizes a given functional J [37j . We consider a functional with a Lagrangian 
depending on the independent variable x, function y and its left Caputo fractional derivative of 
order < a < 1, 



Av] 



L (x, y, c a D a x y) dx = T( a ) a J a b [{b - x^L (x, y, ^D° x 



(15) 



where a J b is seen as the left Riemann-Liouville fractional integral at x = b, subject to the 
boundary conditions y{a) = y a and y(b) = y^, y ai yb S R. For a = 1, J is a functional of the 
classical calculus of variations |44j . In order to obtain a necessary condition for the extremum of 
(Ti"5")) . we use the following two lemmas. 

Lemma 3.1. Let f £ L\([a, b]) and < a < 1. // there is a number e G]a, b] such that 

\f(x)\<c(x-af 
for all x € [a, e] with c > and /3 > —a, then 



Proof. From @ one has 



-du 



< c 



lim a J a x f(x) = 0. 

x— ya+ 



l a {x-u) 1 ~ a ^ T(a + ^ + l) 
for all x 6 [a,e]. Taking into account that a + (3 > 0, it follows that 

fin) 



lim , 

x^a+ J a (X — U) 



1-a 



All 



~ T(a + f3 + 1) a:-»-o+ V ^ 



(16) 



□ 
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Remark 3.1. If < f3 < 1, i/ien function f of Lemma \3.1\ is a Holder function of order fj. If 
j3 = 0, then f is a simple bounded function. 

Lemma 3.2 (The DuBois-Reymond fundamental lemma of the fractional calculus of variations). 
Let g be a differentiable function in [a,b] with g{a) — g{b) — 0, and let f £ Li([a,b}) be such that 
there is a number e £]a, b] with |/(a;)| < c{x — a)* 9 for all x £ [a, e], where c > and (3 > —a are 
constants. Then, 

a J a h (f{xf a lf x g{x)) =0=>f = K, 

where K is a constant. 

Proof. For any constant K we have 

a J a b ({f{x) K) c a D a x g{x)) = a r h (f{x) c a Dlg{ X )) - K a J^D a x g{x) 

= -K a jf a D a x g{x) (17) 
= -K (g(b) - g(a)) 
= 0, 

where we used the fundamental theorem of Caputo calculus and the hypothesis g(a) = g(b). 
Let us choose 

g{x):= a J a x {f{x)-K). (18) 

It can be seen that g(x), defined by (Tl8|) , is a differentiable function in [a, b). Furthermore, from 
(Unj) of Lemma I3TT1 we have g(a) = and, by choosing 

K - 1 T a f - F ( Q + 1 ) ra f 

~ aJt^ (b-a)« aJbh 
we also have g{b) = 0. By inserting (fT8|) into ([T7| . and using ([T0|) . we get 

j_ [ b {_m_Kf 

T(a)J a {b-uf 

Since b — x > for all x e [a, b[, we conclude that f(x) = K for all x € [a, b]. □ 

Remark 3.2. The fractional DuBois-Reymond lemma can be trivially formulated for Riemann- 
Liouville derivatives instead of Caputo ones. Indeed, ^D x g(x) = a D x g(x) because g(a) = 0, as 
commented in Remark \2.4\ Consequently, all results of the next section can also be formulated for 
functionals depending on Riemann-Liouville derivatives. 

4 Fractional Euler— Lagrange Equations 

The next Theorem gives a necessary condition for a function y to be an extremizer of the fractional 
variational functional defined by (|15p . 



a J a h (f(x) - K) c a D a x g{x) = a r h {f{x) ~Kf = — I ^ > > du = 0. 



Theorem 4.1 (The fractional Euler-Lagrange equation in integral form). Let J be a functional 
of the form 

J[y} = J L( K x,y, c a Dly)dx^T{a) a J'i[{b-x) 1 - a L(x,y, c a D a x y) 

defined in the class of functions y £ C ([a, b]) satisfying given boundary conditions y{a) = y a and 
y(b) — yb, and where L £ C 1 ([a, 6[xM 2 ) is differentiable with respect to all of its arguments. If y 
is an extremizer of J, then y satisfies the following fractional Euler-Lagrange integral equation: 



dL(x,yC a D a x y) dL{x,y^D x y) K 

dy + d{CD a x y) (b-xy-° 

for all x £ [a, b[, where K is a constant 



G 



Proof. Let y* give an extremum to (|15p . We define a family of functions 

y(x) = y*(x) + er,(x), (20) 

where e is a constant and r\ € C 1 ([a, b]) is an arbitrary continuously differentiable function satis- 
fying the boundary conditions 77(a) = 77(6) = (weak variations). From (1201) and the boundary 
conditions 77(a) = 77(6) = and y*(a) — y a , y*{b) — y b , it follows that function y is admissible: 
y G C 1 ([a, b]) with y(a) = y a and y(b) = y b - Let the Lagrangian L be C 1 ([a, 6[xR 2 ). Because y* 
is an extremizer of functional J, the Gateaux derivative <5J[y*] needs to be identically null. For 
the functional (fT5|) . 

5J[y*] = lim^ (^j\(x,yC a D a x y) dx - j\{x,y* °D a x y*)dx\ 

' b / nfxi aL(a? ' y *'" g ' g * ) + ^%f ^ L( ^*'°^ V) > l & (21) 
= 0. 

By using ([3]) and relations (0) and (IT5|) . we get 



(22) 



Inserting (|22|) into (l2~Tj) . and using the definition ([2} of left Riemann-Liouville fractional integral, 
we obtain for the first variation the following expression: 

6J\V*] t C D%(x\ ( j ^M^lv*) , dL(x,y-°D a x y-) \ 

= r(a)„ J? 



f v g^jAjg^r) ^(,,7/,^*) X _ Q 



0. 

The fractional Euler-Lagrange equation (fH?|) follows from Lemma 13.21 Note that the hypothesis 
|/(a;)| < \c\(x-a)P with /3 > -a is satisfied because L G C 1 ([a, o[xR 2 ) and y G C 1 ([a, 6]). Indeed, 

in our case /(a;) = ^ J° dL ^ X ' V dv a , D * V } + (6 - a;) 1- " tends to a number c when 

x — >• a. Therefore, /3 = > —a and we are in conditions to apply Lemma 13.21 □ 
It is important to mention that in all previous approaches to the fractional calculus of variations, 
one eliminates the fractional derivative of function 77 appearing in the second integrand of (|21l) 
by performing an integration by parts |37| . However, from the integration by parts (j 13[) . this 
procedure gives an Euler-Lagrange fractional differential equation involving Riemann-Liouville 
derivatives (see [38| and references therein): 

dL{x,y°D a x y) dL(x,yCD a x y) 

dy +x " d{°D a x v) "°- (23) 

In contrast, in the proof of Theorem 14 . 1 1 we take a different procedure by computing an integration 
by parts on the fist integrand instead of the second one. The new fractional DuBois-Reymond 
lemma takes then a prominent role. 

One can obtain the Euler-Lagrange equation ([23| from our result, by taking the right Riemann- 
Liouville derivative to both sides of (fP9l): 



j- r , ja dL ( x ,y>a D xy) 1 r.a dLjx^D^y) _ 1 

dy d(CD x y) (b 
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The Euler-Lagrange equation (|2Uf follows from (|2"41 by Theorem 12. 21 and the well-known equality 
x Dfr(b — x)^ 1 = (see, e.g., Property 2.1 of [3]). A more important consequence from our opti- 
mality condition (|19p . however, is that we can obtain a new fractional Euler-Lagrange differential 
equation involving only Caputo derivatives. 

Theorem 4.2 (The Euler-Lagrange equation with only Caputo derivatives). Consider the problem 
of extremizing (|15[) with a Lagrangian L G C 2 ([a, b] x M 2 ) subject to boundary conditions y(a) = y a 
and y(b) — yt,. If y G C 1 ([a, 6]) is a solution to this problem, then y satisfies the fractional Euler- 
Lagrange differential equation 

dL(x,yCD a x y) c dLjx^D^y) _ 

dy +x 6 d(SDly) - a (25) 

Proof. The optimality condition (|25p is obtained taking the right Caputo fractional derivative 
(|7|) to both sides of (fl~9"l) . Taking the limit x — > b on the left-hand side of (|19p . and having in 

mind that L G C 2 ([a, 6] x M 2 ) and y G C^fa,^), we obtain the quantity ^^fp^j^^ . On 
the other hand, the right-hand side of (TT9l diverges when x — > b if K ^ 0. We conclude that, 
for i G C 2 ([a, 6] x R 2 ), one has if = and — ^g^f^ = (this is a necessary condition 
for the Lagrangian L and the solution y to be nonsingular at x = b). Finally, we obtain (f2"5|) by 
computing the Caputo derivative of the right-hand side of (fT9")) and using Theorem 12.21 □ 
From the well-known relations between the Riemann-Liouville and Caputo fractional deriva- 
tives (see Remark 12. 4[) , the Euler-Lagrange equation in terms of Caputo fractional derivative 
(|25p can be derived from the Euler-Lagrange equation in terms of Riemann-Liouville fractional 



derivative (I23[) by assuming that — ^-^SrTf^ = 0. Here, we do not assume a priori that 



dL{b,y(b),^D^ y(b)) _ s j 10wm g ft as a consequence of the new fractional Euler-Lagrange equation 

in integral form (I19p . The Euler-Lagrange equation (|25p . involving only Caputo derivatives, is 
valid for regular boundary conditions. Consequently, equation (f2"5j) should be better suited to 
applications in physics, sciences and engineering than the fractional Euler-Lagrange equation (|23l) 
with mixed Riemann-Liouville and Caputo fractional derivatives [5]. 



5 Conclusion 

We generalized one of the most important lemmas of the calculus of variations, the DuBois- 
Reymond fundamental lemma of variational calculus, to functionals depending on fractional deriva- 
tives (Lemma I3.2j) . The new lemma enabled us to prove a fractional Euler-Lagrange equation in 
integral form containing fractional derivatives of only one type (Theorem 14. ip . Furthermore, we 
also showed that, when the Lagrangian is a C 2 function, one can then obtain a fractional Euler- 
Lagrange differential equation depending only on Caputo derivatives (Theorem I4.2[) . This is an 
important result because differential equations involving only Caputo derivatives are, in general, 
better suited to applications in physics, sciences and engineering, than a fractional differential 
equation involving both Caputo and Riemann-Liouville derivatives [B]. 
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